Monte Carlo Determination of the Low-Energy Constants of a Spin 1/2 Heisenberg 

Model with Spatial Anisotropy 
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Motivated by the possible mechanism for the pinning of the electronic liquid crystal direction in 
YBCO as proposed in [J], we use the first principles Monte Carlo method to study the spin 1/2 
Heisenberg model with antiferromagnetic couplings Ji and J2 on the square lattice. In particular, 
the low-energy constants spin stiffness ps, staggered magnitization Ms and spin wave velocity c 
are determined by fitting the Monte Carlo data to the predictions of magnon chiral perturbation 
theory. Further, the spin stiffnesses psi and ps2 as a function of the ratio J2/J1 of the couplings are 
investigated in detail. Although we find a good agreement between our results with those obtained 
by the series expansion method in the weakly anisotropic regime, for strong anisotropy we observe 
discrepancies. 
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Introduction. — Understanding the mechanism re- 
sponsible for high-temperature superconductivity in 
cuprate materials remains one of the most active research 
fields in condensed matter physics. Unfortunately, the 
theoretical understanding of the high-Tc materials using 
analytic methods as well as first principles Monte Carlo 
simulations is hindered by the strong electron correla- 
tions in these materials. Despite this difhculty, much 
effort has been devoted to investigating the properties 
of the relevant t-J-type models for the high-Tc cuprates 
0,1 id. Although a conclusive agreement regarding 
the mechanism responsible for the high-Tc phenomena 
has not been reached yet, it is known that the high-Tc 
cuprate superconductors are obtained by doping the an- 
tiferromagnetic insulators with charge carriers. This has 
triggered vigorous studies of undoped and lightly doped 
antiferromagnets. Today, the undoped antiferromagnets 
on the square lattice such as La2Cu04 are among the 
quantitatively best understood condensed matter sys- 
tems. 

Spatially anisotropic Heisenberg models have been 
studied intensely due to their phenomenological impor- 
tance as well as from the perspective of theoretical inter- 
est For example, numerical evidence indicates 
that the anisotropic Heisenberg model with staggered ar- 
rangement of the antiferromagnetic couplings may be- 
long to a new universality class, in contradiction to the 
0(3) universality predictions [lO|. Further, it is argued 
that the Heisenberg model with spatially anisotropic cou- 
plings Ji and J2 is relevant to the newly discovered pin- 
ning effects of the electronic liquid crystal in the under- 
doped cuprate superconductor YBa2Cu306. 45 It 
is observed that the YBa2Cu3 06.45 compound has a tiny 
in-plane lattice anisotropy which is strong enough to pin 
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FIG. 1; The anisotropic Heisenberg model investigated in this 
study. Ji and J2 are the antiferromagnetic couplings in the 
1- and 2-directions, respectively. 



the orientation of the electronic liquid crystal in a partic- 
ular direction. The authors of [l[ demonstrated that the 
in-plane anisotropy of the spin stiffness of the Heisen- 
berg model with spatially anisotropic couplings Ji and 
J2 can provide a possible mechanism for the pinning of 
the electronic liquid crystal direction in YBa2Cu3 06.45. 

Since the anisotropy of the spin stiffness in the spin 
1/2 Heisenberg model with different antiferromagnetic 
couplings Ji and J2 has not been studied in detail be- 
fore with first principles Monte Carlo methods, in this 
letter we perform a Monte Carlo calculation to deter- 
mine the low-energy constants, namely the spin stiff- 
nesses psi and Ps2, staggered magnitization Ms and spin 
wave velocity c. In particular, we investigate the J2/J1- 
dependence of 0=1 and ps2, and find good agreement with 
earlier studies [l| using series expansion methods in the 
weakly anisotropic regime. Our this finding would lead to 
very strong pinning energy per Cu site in YBa2Cu306.45 
as claimed in [l|. However, deviations appear as one 
moves toward strong anisotropy. We argue that the de- 
viations observed between our results and the naive ex- 
pectation might indicate an unexpected behavior of the 
spin-stiffness ps at extremely strong anisotropy. 

Microscopic Models and Corresponding 
Observables. — The Heisenberg model we consider in 
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this study is defined by the Hamilton operator 

X 

where 1 and 2 refer to the two spatial unit- vectors. Fur- 
ther, Ji and J2 in eq. ([T]) are the antiferromagnetic cou- 
plings in the 1- and 2-directions respectively. A physical 
quantity of central interest is the staggered susceptibility 
(corresponding to the third component of the staggered 
magnetization M^) which is given by 

X. = tV ^Tr[Afi'(0)Af3(i)exp(-/3i?)].(2) 
^1^2 Jo ^ 

Here (3 is the inverse temperature, Li and L2 are the 
spatial box sizes in the 1- and 2-direction, respectively, 
and Z = Ti exjp{—(3H) is the partition function. The 
staggered magnetization order parameter Mg is defined 
as Ms = J2xi~^)^^~^^^ Another relevant quantity is 
the uniform susceptibility which is given by 

1 1 

= -TT- / dt -TT[M^{Q)M'{t)eM-m.{3) 
^1^2 Jo ^ 

Here M ~ Sx is the uniform magnetization. Both Xs 
and Xu can be measured very efficiently with the loop- 
cluster algorithm using improved estimators . In par- 
ticular, in the multi-cluster version of the algorithm the 
staggered susceptibility is given in terms of the cluster 
sizes \C\ (which have the dimension of time), i.e. Xs = 
i3l\l2 C^c I^P)- Similarly, the uniform susceptibility 
LTX; (W?) = TTT, (Ec Wt{C)') is given in terms 
of the temporal winding number Wt = Wt{C) which 
is the sum of winding numbers Wt{C) of the loop-clusters 
C around the Euclidean time direction. Similarly, the 
spatial winding numbers are defined by Wi = Wi(C) 
with i £ {1,2}. 

Low-Energy Effective Theory for Magnons. — 
Due to the spontaneous breaking of the SU{2)s spin 
symmetry down to its U{l)s subgroup, the low- 
energy physics of antiferromagnets is governed by two 
massless Goldstone bosons, the antiferromagnetic spin 
waves or magnons. The description of the low-energy 
magnon physics by an effective theory was pioneered by 
Chakravarty, Halperin, and Nelson in [31 ■ A system- 
atic low-energy effective field theory for magnons was 
further developed in [l^, [l^, [l^. The staggered mag- 
netization of an antiferromagnet is described by a unit- 
vector field e{x) in the coset space SU{2)s/U{l)s ~ S'^, 
i.e. e{x) = (ei(x), 62(0;), e3(x)) with e(x)^ = 1. Here 
X = {xi,X2,t) denotes a point in (2+l)-dimensional 
space-time. To leading order, the Euclidean magnon low- 
energy effective action takes the form 

S[e] = / dxi / dx2 / dt (^-^die- die 
Jo Jo Jo ^2 

+ ^d2e ■ d2e + ^dte ■ dte) , (4) 



where the index i € {1,2} labels the two spatial di- 
rections and t refers to the Euclidean time-direction. 
The parameters ps = ^/PsiPs2, Psi and ps2 are the 
spin stiffness in the temporal and spatial directions, re- 
spectively, and c is the spin wave velocity. Rescaling 
x'l = {ps2/Psi)^^*xi and x'2 = {psi/ Ps2)^^'^X2, eq. (g]) can 
be rewritten as 



S[e] = l\x[ l^x', l^t ^(dle- die 

+ late-ate). (5) 



Additionally requiring L[ = L'2 = L we obey the con- 
dition of square area. Notice the effective field theories 
described by eqs. (|4|) and ([5]) are valid as long as the con- 
ditions LiPpsi ^ 1 and Li/3ps2 ^ 1 for i € {1,2} hold, 
which is indeed the case for the set up of this study. Once 
these conditions are satisfied, the low-energy physics of 
the underlying microscopic model can be captured quan- 
titatively by the effective field theory as demonstrated 
in [1^. Further, in the so-called cubical regime (to be 
defined later) which is relevant to our study, the cut-off 
effects appear in the free energy density only at next-to- 
next-to-next-to-leading order (NNNLO). The finite cut- 
off leads to higher-order terms in the effective Lagrangian 
due to the breaking of some symmetries and it introduces 
the cut-off dependence in the Fourier integrals (sums). 
By employing similar arguments as those presented in 
[1§ |. one can show that higher-order corrections to eq. 
(HI) contain four derivatives and the leading cut-off ef- 
fect in the Fourier integrals (sums) enters the free energy 
density only at NNNLO. Therefore eq. ([5]) is sufficient to 
derive up to next-to-next-to- leading order (NNLO) con- 
tributions to the observables considered here. We have 
further verified that the inclusion of NNNLO contribu- 
tions to the relevant observables considered here lead 
to statistically consistent results with those not taking 
such corrections into account. Hence the volume- and 
temperature- dependence of Xs and Xu up to NNLO (to 
be presented below) are sufficient to describe our numer- 
ical data quantitatively, and the finite cut-off effects are 
negligible. Using the above Euclidean action (O, de- 
tailed calculations of a variety of physical quantities in- 
cluding the NNLO contributions have been carried out 
in [1^. Here we only quote the results that are rel- 
evant to our study, namely the finite-temperature and 
finite-volume effects of the staggered susceptibility and 
the uniform susceptibility. The aspect ratio of a spatially 
quadratic space-time box with box size L is characterized 
by / = {Pc/Ly/^ , with which one distinguishes cubical 
space-time volumes with (3c ~ L from cylindrical ones 
with /3c 3> L. In the cubical regime, the volume- and 
temperature-dependence of the staggered susceptibility 
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FIG. 2: Comparison between our numerical results (data 
points) and the theoretical predictions (solid lines) which are 
obtained by using the low-energy parameters from the fits. 
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where Als is the staggered magnetization density. Fi- 
nally the uniform susceptibility takes the form 
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In ^ and ([7]), the functions f3i{l), and which 

only depend on Z, are shape coefficients of the space-time 
box defined in |18| . 

Determination of the Low-Energy Parame- 
ters and Discussions. — In order to determine the 
low-energy constants for the anisotropic Heisenberg 
model given in ([T]), we have performed simulations 
within the range 0.05 < J2/J1 < 1-0. The cubical 
regime is determined by the condition {J2c ^li^)"^ ) ~ 
(Ec^^iC?) - (EcWtiC)^) (which implies (3c ^L). 




FIG. 3: The J2/J1 -dependence of the spin stiffness p^i and 
of the anisotropic Heisenberg model. While the solid circles 
(black) and squares (red) are the Monte Carlo results of psi 
and ps2, respectively^the up and down triangles are the series 
expansion results of [l[ for psi and ps2, respectively. The solid 
lines are added to guide the eye. 



Notice that since J2 < Ji in our simulations, one must 
increase the lattice size Li in order to fulfill the condition 
(Ec Wi (C)2 ) = ( W2iCy ) because eqs. (© and ^ 
are obtained for a (2 -I- l)-dimensional box with equal 
extent in the two spatial directions. Therefore, an inter- 
polation of the data points is required in order to be able 
to use eqs. ^ and ([71). Further, the low-energy parame- 
ters are extracted by fitting the Monte Carlo data to the 
effective field theory predictions. The quality of these fits 
is good as can be seen from figure [5] (the x^/d.o.f for all 
the fits is less than 1.25). Figure [3] shows psi and ps2, 
obtained from the fits, as functions of the ratio of the an- 
tiferromagnetic couplings, J2/J1. The values of psi {Ps2) 
obtained here agree quantitatively with those obtained 
using the series expansion in [l| at J2/J1 = 0.8 and 0.6 
(0.8, 0.6, 0.4, and 0.2). At J2/J1 = 0.4, the value we 
obtained for p^i is only slightly below the corresponding 
series expansion result in [l|. However, sizable devia- 
tions begin to show up for stronger anisotropics. Fur- 
ther, we have not observed the saturation of psi to a 1-D 
limit, namely 0.25Ji as suggested in [l[, even at J2/J1 
as small as 0.05. In particular, psi decreases slightly as 
one moves from J2/J1 = 0.1 to J2/J1 = 0.05, although 
they still agree within statistical errors. Of course, one 
cannot rule out that the anisotropics in J2/J1 considered 
here are still too far away from the regime where this par- 
ticular Heisenberg model can be effectively described by 
its 1-D limit. On the other hand, the Heisenberg model 
considered here and its 1-D limit are two completely dif- 
ferent systems, because spontaneous symmetry breaking 
appears only in 2-D, still ^ = 00 in both cases. Further, 
the low-temperature behavior of Xu in the 1-D system is 
known to be completely different from that of the 2-D 
system [ll, [l^ . Although intuitively one might expect a 
continuous transition of psi, one cannot rule out an un- 
expected behavior of psi as one moves from this Heisen- 
berg model toward its 1-D limit. In particular, since ear- 
lier studies indicate that long-range order already sets in 
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FIG. 4: The J2/Ji-dependence of the spin wave velocity c 
(left) and the staggered magnitization density A4s (right) of 
the anisotropic Heisenberg model. The solid lines are added 
to guide the eye. 

even for infinitesimal small J2/J1 0j[13, HHj it would be 
interesting to consider even stronger anisotropics J2/J1 
than those used in this study to see how p^i approaches 
its 1-D limit. In addition to psi and ps2i we have ob- 
tained M.S and c as functions of J2/J1 as well from the 
fits (figure 111). The values we obtained for Ms agree 
with earlier results in 0], but have much smaller errors 
at strong anisotropics. 

Next, wc would like to turn to discussing the relevance 
of our results to the pinning effect observed empirically in 
YBa2Cu306.45. In [l[ it is argued that the J2/J1 depen- 
dence of the spin-stiffnesses in the spatially anisotropic 
Heisenberg model studied in this work would lead to a 
very strong pinning energy per Cu site (one order of mag- 
nitude larger compared to the corresponding pinning en- 
ergy in La2Cu04). To be more precise, it is the c^uantity 
K which is defined by Ps2l Ps\ = 1 + i^{J2l Ji — 1) hi the 



weak anisotropy regime that results in the claim made in 
Since the spin-stiffnesses calculated here agree with 
those obtained by series expansion in the weak anisotropy 
regime, which in turn implies that our k agrees with that 
in [l| , we conclude that the pinning energy per Cu site is 
indeed very strong. Hence the in-plane anisotropy of the 
spin stiffness of the Heisenberg model with anisotropic 
couplings Ji and J2 can indeed provide a possible mech- 
anism for the pinning of the electronic liquid crystal di- 
rection in YBa2Cu3 06.45. 

Conclusions. — In this note, we have numerically 
studied the Heisenberg model with anisotropic couplings 
Ji and J2 using a loop cluster algorithm. The corespond- 
ing low-energy constants are determined with high pre- 
cision. Further, the J2/'^i-dependence of p^i and Ps2 is 
investigated in detail and our results agree quantitatively 
with those obtained by scries expansion [l[ in the weakly 
anisotropic regime. On the other hand, we observe dis- 
crepancies between our results and series expansion re- 
sults in the strongly anisotropic regime. However, the 
results of our study still lead to very strong pinning en- 
ergy per Cu site in YBa2Cu306 45 which agrees with the 
claim made by the authors in Finally we find that 
an unexpected behavior of psi might be observed as one 
approaches much stronger anisotropy regime than those 
considered in this study. 
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